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Hubble parameter measurement constraints on dark energy 

Omer Farooq 1 , Data Mania 1,2 , and Bharat Ratra 1 

ABSTRACT 



We use 21 Hubble paramet e r versus red s hift d ata p oints, from ISimon et al. 



fl2005f ). lGaztanaga et al.l fl2009[ ) . IStern et al.l (12010 ). and Moresco et all fl2012al ). 
to place constraints on model parameters of cons tant and time- e volving dark 
energy cosmologies. The inclusion of the 8 new iMoresco et al.l (I2012al ) mea- 
surements results in H(z) constraints more restrictive than those derived by 
Chen fc Ratral (j2011bl ). These constraints are now almost as restrictive as those 
that follow fr om current Type l a supernova (SNIa) apparent magnitude versus 
redshift data (jSuzuki et al.ll2012l ). which now more carefully account for system- 
atic uncertainties. This is a remarkable result. We emphasize however that SNIa 
data have been studied for a longer time than the H(z) data, possibly resulting in 
a better estimate of potential systematic errors in the SNIa case. A joint analysis 
of the H(z), baryon acoustic oscillation peak length scale, and SNIa data favors 
a spatially-flat cosmological model currently dominated by a time-independent 
cosmological constant but does not exclude slowly-evolving dark energy. 



1. Introduction 

The expansion rate of the Universe changes with time, initially slowing when matter 
dominated, because of the mutual gravitational attraction of all the matter in it, and more 
recently accelerating. A number of cosmological observations now strongly support the idea 
that the Universe is spatially flat (provided the dark energy density is close to or time 
independent) and is currently undergoing accelerated cosmological expansion. A majority of 
cosmologists consider dark energy to be the cause of this observed accelerated cosmological 
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expansion[j] This dark energy, most simply thought of as a negative pressure subst ance, 
domin a tes the current cosmolog i cal ene r gy budget. For rey iews of dark energy see [Bass 
(120 111 ). Ijimenezl ( 1201 ll ). iLi et al.l (J2011aJ), iBolotin et al.l ( 1201 ll ). and references therein. 



Three observational techniques provide the strongest evi dence for dark energy: SNIa ap 



parent magnitude measurements as a function of redshift (e.g. , ISullivan et al.ll2011t ISuzuki et al. 



2012 



i et allboilrJ : iBarreira fc Avelinoll2011l); c osmic microwave background (CMB) anisotropy 



data (IPodariu et al.ll2001bl ; iKomatsu et al.ll201ll . a nd references therein ) combined with low 



estimates of the cosmological mass density (e.g., IChen fc Ratral 120031 ) . provided the dark 
energy density is close to or tim e independent; and baryon acoustic oscillation (BAO) pe ak 
length scale measurements (e.g., iBeutler et al.l 1201 It iBlake et al.l 1201 ll ; iMehta et al.ll2012l ). 



The "standard" model of cosmology is the spatially- flat ACDM model (|Peebleslll984l ). In 
this model about 73% of the current energy budget is dark energy, Einstein's cosmological 
constant A. Non-relativistic cold dark matter (CDM) is the next largest contributer to 
the energy budget (around 23%), followed by non-r elativistic baryonic mat ter (about 5%). 
For reviews of the standard model of cosmology see iRatra fc Vogeleyl ( 120081 ) and references 
therein. It has been known for a while that the ACDM model is reasonably consistent 



with most observ ations (see, e.g., IJassal et al.l I201CH : IWilson et al.l 120061 ; iDavis et al.l 12007 



Allen et al. 2008 . for early indications) |§ In the ACDM model the dark energy density is 



constant in time and does not vary in space. 

Although most predictions of the ACDM model are reasonably consistent with the 
measurements, the ACDM model has some curious features. For instance, the measured 
cosmological constant energy density is 120 orders of magnitude smaller than the energy 
density naively expected from quantum field theory (this is known as the fine-tuning puzzle). 
A second curiosity is what is known as the coincidence puzzle: the energy density of a 
cosmological constant, p^, is independent of time, but that of matter, p m , decreases with 
time during the cosmological expansion, so it is curious why we (observers) happen to live at 
this (apparently) special time, when the dark energy and the non-relativistic matter energy 
densities are of comparable magnitude. 



1 Some cosmologists instead view these observations as an indication that general relativity needs t o 
be modified on the s e large length scales. For r e cent re views of modifie d gravity see iTsuiikawal (|2010n . 



Bolotin et al.l (|201ll ). ICapozziello fc De Laurentisl (|201lh . IStarkmanl (|201ll ). and references therein. In this 
paper we assume that general relativity provides an adequate description of gravitation on cosmological 
scales. 



2 Note, however 

model, might need modification (see Peebles fc Ratra 2003 ; Perivolaropoulosl 2010 . and references therein). 



that the "standard" CDM structure formation model^which is assumed in the ACDM 
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These puzzles co uld be partially resolve d if the dark energy den sity is a slowly decreas- 
ing function of time (IPeebles fc Ratralll988l ; iRatra &: Peebles! Il988l ). In this case the dark 
energy density will remain comparable to the non-relativistic matter density for a longe r 
time. For recen t discu s sions of time-varying dark energy models, s ee iBauer et al.l (120111). 



Basilakos ct al. 



Cai et al.l (120121 ) 



Chimento et al.l (120111). iGranda et al.l (120111). iGarcia-Bellido et all ( 120111) . 
(l2012[ ). IShevkhi et al.l (boAlBrax fc Davisl d2012h . lHollenstein et al.l EoV$ \. 
and references therein. In this paper we will consider two dark energy models (with dark 
energy being either a cosmological constant or a slowly-evolving scalar field 0) as well as a 
dark energy parameterization. 

In the ACDM model, time-independent dark energy density (the cosmological constant 
A) is modeled as a spatially homogeneous fluid with equation of state p\ = —p\. Here 
and pa are the fluid pressure and energy density. In describing slowly- decreasing dark 
energy density much use has been made of a parameterization known as XCDM. Here dark 
energy is modeled as a spatially homogeneous X-fluid with equation of state px = wxPx- 
The equation of state parameter wx < — 1/3 is independent of time and px and px are the 
pressure and energy density of the X-fluid. When wx = — 1 the XCDM parameterization 
reduces to the complete and consistent ACDM model. For any other value of wx < — 1/3 
the XCDM parameterization is incomple te as it cannot describe spatial inhomogeneities (see, 
e.g. lRatralll991t iPodariu &: Ratral l2000l ) . For computational simplicity, in the XCDM case 



we assume a spatially-flat cosmological model. 

The 0CDM mod el is the simplest, consi s tent and complete mod el of slowly- decreasing 
dark energy density ( IPeebles fc Ratral Il988l ; IRatra fc Peebles! Il988l ). In this model dark 
energy is modeled as a scalar field, 0, with a gradually decreasing (in 0) potential energy 
density V((p). Here we assume an i nverse power-law pote ntial energy density V(4>) oc <p~ a , 
where a is a nonnegative constant ( IPeebles fc Ratral Il988l ). When a = the 0CDM model 
reduces to the corresponding ACDM case. For computational simplicity, we again only 
consider the spatially-flat cosmological case for 0CDM. 

As mentioned above, for some time now, most observational constraints have been rea- 
sonably consistent with the predictions of the "standard" spatially-flat ACDM model. CMB 
anisotropy, SNIa, and BAO measurements provide the strongest support for this conclusion. 
However, the error bars associated with these three types of data are still too large to allow 
for a significant observational discrimination between the ACDM model and the two simple 
time-varying dark energy models discussed above. This is one motivation for considering 
additional kinds of data. 



If the constraints from the new data differ considerably from the old ones, this could 
mean that at least one of the data sets had an undetected systematic error, or it could mean 
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that the model being tested is observationally inconsistent. Either of these is an important 
result. On the other hand, if the constraints from the new and the old data are consistent, 
then a joint analysis of all the data could result in tighter constraints, and so might result 
in significantly discriminating between constant and time-varying dark energy models. 



Other measurements that have been used to constrain cosmologic al parameter^ in 



elude galaxy cluster ga s mass fraction as a function of redshift (e.g.. lAllen et al.l 12008 



Samushia & Ratrall2008l : Ettori et al.ll2009l ; iTong & Noh 



ter and other large-sca l e structure properties (ICampanelli et al.l 120121 ; 



2011 



Lu et a 



20111). galaxy clus 



De Boni et al. 2011 



Mortonson et al.l l201ll ; bevi et al.l l201ll ; IWaneJ I2OI2L and references therein), gamma-ray 
burst luminosity d i stanc e as a function of redshift (e.g., Samushia &: RatralboiO ; Wang &: Dai 



2011 



iBusti et al 



Dantas et al 



2011 



2012 ), lookback tim e as a function of redshift ( jSamushia et al 



2010 



Tonoiu et al.l 1201 ll . and refe rences therein), HI 



paren t magnitude as a function of redshift (e.g., P lionis et al. 



2012 1), angular size a s a function of redshift (e.g., iGuerra et al 



Chen fc Ratra 120121) , and strong gray itational lensing (jChae et al.l 12004 ; iLee fc Ng 



2010 



starburst galaxy ap- 



2011 



Mania fc Ratra 



20001; Bonamente et al. 



2006 



2007 



Biesiada et al. 2010 ; Zhang fc Wu 2010 . and references therein) |j| Of particula r interest to 



us 



here are measurements o 



2003 



Samushia fc Ratra 



the Hubble parameter as a f unction of redshift (e. g^ 



2006; Samushia et al. 2007; Sen fc Scherrer 



Kumar Il20l3 " IWang fc Zhangl 



2008 



Jimenez et al. 



Pan et al 



201 



0; 



Chen fc Ratra 



2012 



20 



Seikel et al. 



lb 



2011 



Duan et al. 



2011 



Bilicki et al. 



20121 ). While the constraints from these data are typically less restrictive 
than those derived from the SNIa, CMB anisotropy, and BAO data, both types of measure- 
ments result in largely compatible constraints that generally support a currently accelerating 
cosmological expansion. This provides confidence that the broad outlines of a "standard" 
cosmological model are now in place. 



In this paper w e use t he 21 H(z) measurements oflSimon et al.l (120051 ) . Gaztanaga et al. 
J2009I ). Istern et all (bold ), and iMoresco et all (boi2ah (listed in Table 1)| to constrain the 



ACDM and (ftCDM m odels and the XCDM parametrization. The inclusion of the 8 new 
Moresco et al.l (j2012al ) measurements (with smaller error bars compared to t he earlier data) 



in the analysis results in tighter constraints than those recently derived by IChen fc Ratra 



3 For reviews see Albrecht et al. ( 20061 ). Weinberg et al. ( 2012 ). and references therein. 

4 Future space-based SNIa and BAO-like meassuremen ts fe.g.. |Podariu et al 1 l2001at I Samushia et al. I l201ll 



Sartoris et al 



20121 



( Jennings et al 



Basse et al 



2011 



2012t 



Pavlov 



van de Wevgaert et al 



et al 



2011 



2012h. as well as m easurements based on new techniques 



Ziaeepour 1 120121 and references therein) should soon 



provide interesting constraints on cosmological parameters. 



5 We do not include the 4 recent IZhang et al.l (|2012T) H (z) measurements as they have somewhat larger 
error bars and do not affect our results. 



- 5 - 



z 


H{z) 




Reference 




(km s" 1 Mpc - 1 ) 


(km s" 1 Mpc - 1 ) 




0.090 


69 


12 


1 


0.170 


83 


8 


1 


0.179 


75 


4 


4 


0.199 


75 


5 


4 


0.240 


79.69 


2.65 


2 


0.270 


77 


14 


1 


0.352 


83 


14 


4 


0.400 


95 


17 


1 


0.430 


86.45 


3.68 


2 


0.480 


97 


62 


3 


0.593 


104 


13 


4 


0.680 


92 


8 


4 


0.781 


105 


12 


4 


0.875 


125 


17 


4 


0.880 


90 


40 


3 


0.900 


117 


23 


1 


1.037 


154 


20 


4 


1.300 


168 


17 


1 


1.430 


177 


18 


1 


1.530 


140 


14 


1 


1.750 


202 


40 


1 



Last column reference numbers: 1. 



Simon et al. 


(2012a 





(12009( 1. 3. IStern et all (l2010f ). 4. Moresco et al. 
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( l2011bl ) from the previous largest set of H(z) measurements considered. The new H(z) data 
constraints derived here are compatible with cosmological parameter constraints determined 
by other techniques. For the first time, these H(z) limit s are almost a s cons training as those 
derived from the most recent SNIa data compilation of ISuzuki et al.l (120121 ) . In addition to 
the tighter H(z) limits resulting from the new data, this is p artially also a conse quence of the 
fact that a more careful analysis of the SNIa measurements (ISuzuki et al.ll2012l ) has resulted 
in a larger systematic error estimate and thus weaker SNIa constraints. We emphasize that 
the study of H(z) data is much less mature than that of SNIa apparent magnitude data, so 
there is the possibility that future H(z) error bars might be larger than what we have used 
in our analysis here. In addition to deriving if (z)-data only constraints, we also use these 
H(z) data in combination with recent BAO and SNIa measurements to jointly constrain 
cosmological parameters in these models j§ Adding the H (z) data tightens the constraints, 
somewhat significantly in some parts of parameter space for some of the models we study. 

Our paper is organized as follows. In Sec. [2] we present the basic equations of the three 
dark energy models we consider. Constraints from the H(z) data are derived in Sec. [3j In 
Sec. H] we determine constraints from recent SNIa apparent magnitude data. In Sec. we 
derive constraints from recent BAO data. Joint constraints on cosmological parameters, from 
a combined analysis of the three data sets, for the three models we consider, are presented 
in Sec. [6j We conclude in Sec. 



2. Dark energy models 

In this section we summarize properties of the two models (ACDM and 0CDM) and the 
one parametrization (XCDM) we use in our analyses of the data. 

To determine how the Hubble parameter H(z) evolves in these models, we start from 
the Einstein equation of general relativity 

Rfiu - -g^R = %kGT^ u - Ag^ lu . (1) 

Here g^ v is the metric tensor, R^ u and R are the Ricci tensor and scalar, T^ u is the energy- 
momentum tensor of any matter present, A is the cosmological constant, and G is the 
Newtonian gravitational constant. 

The energy- momentum tensor for an ideal fluid is T^ v = diag(p,p, p,p), where p is the 



6 See Moresco et al. ( 2012a ) and Wang et al. ( 2012 ) for analyses that use most of these H(z) data in 
conjunction with CMB anisotropy and other data to constrain cosmological parameters. 
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energy density and p the pressure. Assuming spatial homogeneity, the Einstein equation 
reduces to the two independent Friedmann equations 

a\ 2 8nG A K 2 

P+o--^ (2) 



3 3 a 

a 4:iiG A 

a 3 3 

Here a(t) is the cosmological scale factor, an overdot denotes a derivative with respect to 
time, and K 2 represents the curvature of the spatial hypersurfaces. These equations, in 
conjunction with the equation of state, 

p = p (p) = ujp, (4) 

where u is the dimensionless equation-of-state parameter (with cj = — 1 corresponding to a 
cosmological constant and u < —1/3 corresponding to the XCDM parametrization), govern 
the evolution of the scale factor and matter densities. 

Taking the time derivative of Eq. ([2]) and putting it in Eq. (J3J) and then using Eq. (TjJ 
yields the energy conservation equation 

p = -3-(p + p) = -3p-(l + u) (5) 
a a 

For a non-relativistic gas (matter) u = u m = and p m oc a -3 , and for a cosmological constant 
uj = cj\ = — 1 and pa = A/(87rG)= constant (p\=0). Solving Eq. (jHJ), the time- dependent 
energy density is 

where po and ao are the current values of the fluid energy density and the scale factor. If there 
are a number of different species of non-interacting particles, then Eq. OH]) holds separately 
for each of them. 

The ratio a(t)/a(t) = H(t) is called the Hubble parameter. The present value of the 
Hubble parameter is known as the Hubble constant and is denoted by H . Defining the 
redshift z = a /a — 1, and the present value of the density parameters, 

_ SnGpp -K 2 A 

in the ACDM model we can rewrite Eq. §Z§ as 

H 2 {z- H , p) = Hi [Q m0 (l + zf + n A + (1 - n m0 - fi A )(l + zf] , (8) 
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where we have made use of Qko = 1 — H m o- Q\. This is the Friedmann equation of the ACDM 
model with spatial curvature. In this model the cosmological parameters p = (Q m o,Q\). 
Here Q m o is the non-relativistic (baryonic and cold dark) matter energy density parameter 
at the present time. Below we shall have need for the dimensionless Hubble parameter 
E{z) = H(z)/H . 

It has become popular to parametrize time-varying dark energy as a spatially homoge- 
neous X-fiuid, with a constant equation of state parameter w x = Px/px < —1/3. With this 
XCDM parametrization the Friedmann equation takes the form 

E\Z- H , P ) = # 2 [On0(l + Zf + (1 - On<))(l + zf {1+ ^\ (9) 

where for computational simplicity we consider only flat spatial hypersurfaces, and the model 
parameters p = (Q m o,ux)- The XCDM parametrization is incomplete, as it cannot describe 
the evolution of energy density inhomogeneities. 

The simplest complete and consistent dynamical dark energy model is 0CDM. In this 
model dark energy is a slowly-rolling scalar field <fi with an, e.g., inverse-power-law potential 
energy density V(4>) = nm 2 ) (j)~ a where m p = 1/a/G is the Planck mass and a is a non- 
negative free parameter that determines k. The scalar field part of the 0CDM model action 
is 

ml f , / 1 



s = / v 73 ^ ( g^d^dvt - Km 2 p( f)- a ) <fx, (10) 



16-7T J ° \2 

with corresponding scalar field equation of motion 

4> + 3-0 - Kam 2 J- {a+1) = 0. (11) 
a p 

In the spatially-flat case the Friedmann equation is 

H 2 (z; H , p) = ^(p m + p*) = H 2 [n m0 (l + zf + %(z, «)], (12) 

with scalar field energy density given by 

ml ( 1 



I6i{^ 2 + Km >" (13) 
Solving the coupled d ifferential Eqs. ( ITT]) — (fT3|) . with the initial conditions described in 



Peebles fc Ratra 



(119881 ). allows for a numerical computation of the Hubble parameter H(z). 



In this case the model parameter set is p = (Q m o, o). 
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Constraints from the H(z) data 



We use 21 i ndependent H(z) da ta points (jSimon et al.l 120051 ; iGaztanaga et al.l 12009 



Stern et al.ll2010l ; lMoresco et al.ll2012al ). listed in Table 1, to constrain cosmological model pa- 
rameters. The observational data consist of measurements of the Hubble parameter H Q ^ s (zi) 
at redshifts z i7 with the corresponding one standard deviation uncertainties <7j. 

To constrain cosmological parameters p of the models of interest we compute the Xh 
function 



21 



X 2 H (H ,P) = J2 



[H th (zi; H , p) - H oha (zi)Y 



i=l 



07 



(14) 



where H th (zi] H , p) is the model-predicted value of the Hubble parameter. As discussed in 
Sec. El H th (zi; H Q , p) = H E(z; p), so from Eq. (JTU) we find 



21 



X 2 H (Ho,p) = # 2 E 



E 2 (z i]P ) 



21 



i=l 



(J 



H ohs (zi)E(zi;p) 



21 



i=l 



at 



(15) 



i=i 



Xh depends on the model parameters p as well as on the nuisance parameter H whose 
value is not known exactly. We assume that the distribution of i^o is a Gaussian with one 
standard deviation width cr# and mean H . We can then build the posterior likelihood 
function £#(p) that depends only on the p by integrating the prod uct of exp(—Y 2 H /2) and 
the H prior likelihood function exp[— (H — H ) 2 / (2a H )] (see, e.g-. lGanga et al.lll997l ). 



27U7 k o 



e -X 2 H (Ho,p)/2 e -(H -H f/(2aj lo ) dH ^ 



(16) 



Defining 



21 



a 



L> TJ.. 



E 2 (z l]P ) 



'H — i a i "H — i "i "lt„ 

the integral can be expressed in terms of the error function[_] 



H A H ohs (z i )E(z i ;p) H 2 A H 2 ohs ( Zi ) 

( 3 = ^2- + l^ ~2 > ^ = ^2- + l^ 

<J " i=l a i a H i=1 



a? 



(17) 



2ja a 2 Ho 



exp 



2 



1 + erf 







7 erf (x) 
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0.4 0.6 0.8 1.0 



Fig. 1. — Solid lines shows 1, 2, and 3 a constraint contours for the ACDM model from the 
H(z) data. The left panel is for the Ho = 68 ±2.8 km s _1 Mpc" 1 prior and the right panel is 
for the Hq = 73.8 ± 2.4 km s _1 M pc" 1 one. Thin dot-das hed lines in the left panel are 1, 2, 
and 3 a contours reproduced from lChen &: Ratral (l2011bl ). where the prior is Hq = 68 ± 3.5 
km s _1 Mpc -1 ; the empty circle is the corresponding best-fit point. The dashed diagonal lines 
correspond to spatially-flat models, the dotted lines demarcate zero-acceleration models, and 
the shaded area in the upper left-hand corners are the region for which there is no big bang. 
The filled black circles correspond to best-fit points. For quantitative details see Table EJ 



We maximize the likelihood £#(p), or equivalently minimize Xh(p) — ~ 21n£#(p), 
with respect to the parameters p to find the best-fit parameter values p . In the models we 
consider x\ depends on two parameters. We define la, 2a, and 3a confidence intervals as 
two-dimensional parameter sets bounded by Xh(p) = Xh(Po) + 2.3, Xh(p) — Xh(Po) + 6.17, 
and Xh(p) = Xh(Po) + H-8, respectively. 

Even though the precision of measurements of the Hubble constant have greatly im- 
proved over the last decade, the concomitant improvement in the precision of other cosmo- 
logical measurements means that in some cases the Hubble constant uncert ainty still signif- 



i cantly affects cosmological parameter estimation. For a recent example see ICalabrese et al. 



(120121 ). The values of H Q ± au that we use in this paper are 68 ± 2.8 km s 1 Mpc 1 a nd 



73.8 ± 2.4 km s 1 Mpc 1 . The first is from a me dian statistics analysis (IGott et al.ll200ll ) of 



(Ch en fe Ratra 



553 measurements of H 
for over a decade now ( Gott et al. 



2001 



2011ah: this estim ate has been remarkably stable 
The second value is the most 



Chen et al 



precise recent one, based on HST measurements (IRiess et al.ll201ll ). Other recent measure- 



2003h. 
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Fig. 2. — Solid lines shows 1, 2, and 3 a constraint contours for the XCDM parametrization 
from the H(z) data. The left panel is for the Hq = 68 ± 2.8 km s _1 Mpc -1 prior and the 
right panel is for the H = 73.8 ± 2.4 km s -1 Mpc -1 one. Thin dot - dashed lines in the 
left panel are 1, 2, and 3 a contours reproduced from IChen fc Ratral (j2011bl ). where the 
prior is Hq = 68 ± 3.5 km s -1 Mpc -1 ; the empty circle is the corresponding best-fit point. 
The dashed horizontal lines at wx = — 1 correspond to spatially-flat ACDM models and the 
curved dotted lines demarcate zero-acceleration models. The filled black circles correspond 
to best-fit points. For quantitative details see Table |2j 



ments are not incons i stent with at least one of the two values we use as a prior (see, e.g., 



Freedman et al.l 120121 ; ISorce et al.l 12012c iTammann fc Reindll 120121 ) . 



Figures [T]— [3] show the constraints from the H(z) data for the three dark energy models 
we consider, and for the two different Hq priors. Table [2] l ists the best fit parameter values. 
Comparing these plots with Figs. 1 — 3 of IChen fc Ratral (j2011bl ). whose 1, 2 and 3 a con- 
straint contours are reproduced here as dot-dashed lines in the left panels of Figs [Q— El we 
see that the contours derived from the new da ta are more constrainin g, by about a standard 
deviation, because of the 8 new, more precise, iMoresco et al.l ( )2012al ) data points used here. 
On comparing the left and right panels in these three figures, we see that the constraint 
contours are quite sensitive to the value of H used, as well as to the uncertainty associated 
with the Hubble constant measurement. 
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Fig. 3. — Solid lines shows 1, 2, and 3 a constraint contours for the 0CDM model from the 
H(z) data. The left panel is for the Hq = 68 ±2.8 km s _1 Mpc -1 prior and the right panel is 
for the H = 73.8 ± 2.4 km s _1 M pc -1 one. Thin dot-das hed lines in the left panel are 1, 2, 
and 3 a contours reproduced from lChen fc Ratral (l2011bl ). where the prior is Hq = 68 ± 3.5 
km s _1 Mpc -1 ; the empty circle is the corresponding best-fit point. The horizontal axes at 
a = correspond to spatially-flat ACDM models and the curved dotted lines demarcate zero- 
acceleration models. The filled black circles correspond to best-fit points. For quantitative 
details see Table [2j 



Constraints from the SNIa data 



While the H(z) data provide tight constraints on a linear combination of cosmological 
parameters, the very elongated constraint contours of Figs. Q3— [3] imply that these data alone 
cannot significantly discriminate between cosmological models. To tighten the constraints 
we must add other data to the mix. 



The second set of data that we use are the Type la supernova data from the lSuzuki et al. 



(120121 ) Union2.1 compilation of 580 SNIa distance modulus ^ bs( z i) measurements at mea- 
sured redshifts Zi (covering the redshift range of 0.015 to 1.414) with associated one standard 
deviation uncertainties (Tj. The predicted distance modulus is 



/ithte; H , p) = 5 log 10 (3000 y{z){\ + z)) + 25 - 5 \og w (h) : 



(19) 



— 
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H(z) 




SNela 




RAO 


Model and prior 


y 2 ■ 

A mm 


R F P 

J_J . JL . J. 


v 2 ■ 

A mm 


R F P 


2 

Xmin 


R F P 


7 

h = 


ACDM 
0.68 ± 0.028 


14.6 


^m0 — 0.28 


545 


O™o = 0.29 


5.5 


^m0 = 0.27 


h — 


ACDM 
n 7Qs -i- n no/I 

U. ( OO It u.uz^t 


14.6 


fi m0 = 0.42 
Qa — 97 


fi A =0.69 


fi A =0.87 


h = 


XCDM 
0.68 ±0.028 


14.6 


^mO — 0.31 

uiv — —0 Q4 


545 


n m0 = 0.29 


5.5 


fi m o = 0.27 


h = 


XCDM 
0.738 ± 0.024 


14.6 


n m0 = o.3o 

ojx = —1.3 


Wx = -0.99 


uj x = -1.21 


h = 


0CDM 
0.68 ±0.028 


14.6 


fi m0 = 0.30 
a = 0.25 


545 


Q m0 = 0.27 


5.9 


Q m0 = 0.30 


h = 


0CDM 
0.738 ± 0.024 


15.6 


O m0 = 0.27 
a = 0.00 


a = 0.20 


a = 0.00 



Table 2: The minimum value of v 2 and the corresponding best-fit points (R.F.P) which maxi- 
mize the likelihood for the three individual data sets. The SNIa values are for the case includ- 
ing systematic errors. Ignoring SNIa systematic errors, for the ACDM model x|at(Po) = 562, 
at (fi m0 , n A ) = (0.28, 0.73); for the XCDM case xijv(Po) = 562 at (Q m0 , u x ) = (0.28, -1.01); 
and for the 0CDM model vjw(Po) = 562, at (Q m0 ,a) = (0.27,0.05). 



where H = lOO/i km s 1 Mpc 1 and y(z) is the dimensionless coordinate distance, 

f dz' 
J E(z') 



y(z) 



agHu ■ I K 
K blU \a H J E(z') 



K 2 = 



(20) 



sinh ( ^mr I W) ) K2 <°- 



As the SNIa distance modulus measurements /i bs are correlated, v 2 is defined as 

X 2 SN (h,p) = Afi T C- 1 Afi. (21) 

Here A/x is a vector of differences A/Zj = /ith(^; H , p) — ^obs(zi), and C~ x is the inverse of 
the 580 by 580 Union 2.1 compilation covariance matrix. In index notation, 



XSN 



(K P) = \M0 ~ 51 °gl0 /i - A*obs] Q (C \Hq - 5l0g 1D /l - /iobs]^ • 



(22) 
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The covariance matrix is symmetric so this can be written as 

X 2 SN (h,p) = A(p) - 10S(p)log 10 (/i) + 25C[log 10 (/ i )] 2 

where 

Mv) = XX/^O - Hobs)a (C _1 )o/3 (/*0 - yUobs)/? 
a,/3 

5 (P) = ZXA*0 - /iobs)a E( C_1 )a/3 

c = E(c- 1 W 

The corresponding likelihood function, when considering a flat H prior, is 



oo 



(23) 



(24) 



(25) 



Defining 



tf _ 25C ^ _ £(p)lnlO 



2(lnl0) : 



the above integral takes the form 



£-Sn(p) 

The /i-independent 

xIjv(p) = 



cxp 



- 2 I A(P) 



5C 



B 2 (P) „ 1 
C 25* 



-2 ln£ SiV (p) = A(p) 



B 2 (p) 21n(10) 



5(p)-Q, 



C 5C 

where Q is a constant that does not depend on the model parameters p 



625 C 2 



25 C 



(26) 



(27) 



and so can be ignored. We minimize Xsjv(p) with respect to the model parameters p to find 
the best-fit parameter values p and constraint contours. 

Figure H] shows constraints from the SNIa data on the three dark energy models we 
consider here. For the ACDM model and the XCDM parametrizat ion the constraints shown 
in Fig. H] are in very good agreement with those in Figs. 5 and 6 of ISuzuki et all ( 120121 ). The 
</>CDM model SNIa data constraints shown in Fig. H] have not previously been computed. 
Comparing the SNIa constraints of Fig. H] to those which follow from the H(z) data, Figs. 
[33— [U it is clear that SNIa data provide tighter constraints on the ACDM model. For the 




Fig. 4. — Thick solid (dot-dashed) lines are 1, 2, and 3— a constraint contours from SNIa 
data with (without) systematic errors. Filled (open) circles demarcate likelihood maxima 
for the case of data with (without) systematic errors. The top left plot is for the ACDM 
model, the top right plot is for the XCDM parametrization, and the bottom one is for the 
0CDM model. For quantitative details see Table EJ 
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XCDM case both SNIa data and H{z) data provide approximately similar constraints, while 
the SNIa constraints are somewhat more restrictive than the H(z) ones for the 0CDM model. 
However, in general, the SNIa constraints are not very significantly more restrictive than the 
H(z) constraints, which is a remarkable result. It is also reassuring that both data favor 
approximately similar regions of parameters space, for all three models we consider. However, 
given that the degeneracy in parameter space is similar for the H(z) and SNIa data, a joint 
analysis of just these two data sets is unlikely to greatly improve the constraints. 



5. Constraints from the BAO data 



In an attempt to further tighten the cosmological parameter constraints, we now include 
BAO data in the analysis. To const rain cosmological parameters using BAO data we follow 
the procedure of iBlake et al.l (120111 ). To derive the BAO constraints we make use of the 
distance parameter Dy(z), a combination of the angular diameter distance and the Hubble 
parameter, given by 

\l + zfd A {zf CZ 



D v {z) 

Here d A (z) is the angular diameter distance 

d A {z) = 



Hz) 



1/3 



y{z) 



H (l + z) 

where y(z) is the dimensionless coordinate distance given in Eq. (T2"Uj) . 



(2f 



(29) 



We use measurements of the acoustic parameter A(z) fromlBlake et al.l (120111) . where 
the theoretically-predicted A th (z) is given in Eq. (5) of lEisenstein et al.l ( 120051 ). 



Ath(z) 



100 D v (z) y/VL^h 1 



(30) 



Using Eqs. (|25jl— (|30j) we have 



Ath(z) 



y\z) 

z 2 E(z) 



1/3 



(31) 



which is h independent and where E(z) is defined in Sec. [2j 



Using the WiggleZ A Q i, s (z) data from Table 3 of IBlake et al.l (120 111 ), we compute 

x ^(p) = AA r (C A J- 1 AA. (32) 
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Fig. 5.— 1, 2, and 3 a constraint contours from the BAO data. Filled circles denote 
likelihood maxima. The top left plot is for the ACDM model, the top right one is for the 
XCDM parametrization, and the bottom plot is for the 0CDM model. For quantitative 
details see Table El 
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Here AA is a vector consisting of differences AAi = A t h(zj - p) — A^iz i ) and (Ca z ) 1 is the 
inverse of the 3 by 3 covariance matrix given in Table 3 of iBlake et all (1201 if ). 



We also use the 6dF GS and S PSS data, th r ee me asurements from lBeutler et all (1201 if ) 
and |Per rival et al.l (120101 ) , listed in IBlake et all (j201l[ ) . In this case the distilled parameter 



d t h (z) 



r s {zd) 
Dv(z) 



(33) 



where r s (zd) is the sound horizon at the drag epoch, is given in Eq. (6) of lEisenstein Hu 



( 1998 ). The correlation coefficients for this case are also given in Table 3 of IBlake et al. 
(1201 ll ) . Using the covariance matrix we define 



X ^(^,p) = Ad r (C d J- 1 Ad 



(34) 



where Ad is a vect or consisting of di fferences Adi = d^zf, h, p) —d h s (zi) and Cd z is the the 
covariance matrix (IBlake et al.ll201ll ). We then marginalize over a flat prior for Hq to get 



xl (p) 



-2 In 



-X 2 dz (h, P )/2 dh 



(35) 



Since Xa z (p) an d Xd (p) correspond to independent data, the combined BAO data 



XbaoKP) 



X 2 a z 



xL(p) 



(36) 



We can maximize the likelihood by minimizing Xbao(p) with respect to the model param- 
eters p to get best-fit parameter values po and constraint contours. Figure show the 
constraints from the BAO data on the three dark energy models we consider here. The 
XCDM parametriza tion constraints sh own in this figure are in good agreement with those 
shown in Fig. 13 of IBlake et all (120111 ). The constraints shown in the other two panels of 
Fig. have not previously been computed. Comparing to the H(z) and SNIa constraint 
contours of Figs. [I]— [H we see that the BAO contours are also very elongated, although 
largely orthogonal to the H{z) and SNIa ones. Consequently, a joint analysis of these data 
will result in significantly tighter constraints than those derived using any one of these data 
sets. 



6. Joint constraints 

To constrain cosmological parameters from a joint analysis of the H(z), SNIa, and BAO 
data we compute 

X 2 (P) = Xh(p) + Xsn(p) + X 2 bao(p) (37) 
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Fig. 6.— Thick solid (dot-dashed) lines are 1, 2, and 3 cr constraint contours for the 
ACDM model from a joint analysis of the BAO and SNIa (with systematic errors) data, 
with (without) the H(z) data. The full (empty) circle marks the best-fit point determined 
from the joint analysis with (without) the H(z) data. The dotted sloping line corresponds 
to spatially-flat ACDM models. In the left panel we use the Hq = 68 ± 2.8 km s _1 Mpc -1 
prior while the right panel is for the H = 73.8 ± 2.4 km s -1 Mpc^ 1 case. For quantitative 
details see Table [31 



for each of the three cosmo logical models considered here. We minimize x 2 (p) with respect 
to model parameters p to get best-fit parameter values po and constraint contours. 

Figures EHS] show constraints on the cosmo logical parameters for the ACDM and </>CDM 
models and the XCDM parametrization, from a joint analysis of the BAO and SNIa data, as 
well as from a joint analysis of the BAO, SNIa and H(z) data. Table 3 lists information about 
best-fit parameter values. Including the H(z) data in the analysis tightens the constraints 
by more than one standard deviation, in parts of the parameter spaces. 

Adding the H(z) data for the H ± a Ho = 68 ± 2.8 km s" 1 Mpc" 1 prior case improved 
the constraints most significantly in the ACDM case (by more than 1 a on £7a in parts of 
parameter space), Fig. [61 and least significantly for the 0CDM model, Fig. [HI For the case of 
the Hq±<jh = 73.8 ±2.4 km s -1 Mpc -1 prior, adding H(z) again tightens up the constraints 
the most for the ACDM model (by more than 1 a on Fig. [HI and least so for the XCDM 
parametrization, Fig. [71 
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Fig. 7. — Thick solid (dot-dashed) lines are 1, 2, and 3 a constraint contours for the XCDM 
parametrization from a joint analysis of the BAO and SNIa (with systematic errors) data, 
with (without) the H(z) data. The full (empty) circle marks the best-fit point determined 
from the joint analysis with (without) the H(z) data. The dotted horizontal line at wx = — 1 
corresponds to spatially-flat ACDM models. In the left panel we use the Hq = 68 ± 2.8 km 
s _1 Mpc -1 prior while the right panel is for the H = 73.8 ± 2.4 km s _1 Mpc -1 case. For 
quantitative details see Table [3j 

Figures 191— TTTl show the constraints on the cosmological parameters of the three models, 
from a joint analysis of the BAO and H(z) data, as well as from a joint analysis of the three 
data sets. Table 3 lists the best-fit parameter values. Comparing these figures to Figs. [6]— [8] 
allows for a comparison between the discriminating power of the SNIa and H(z) data. 

Figure [9] shows that adding SNIa data to the H(z) and BAO data combination for the 
Hq ± a Ho = 68 ± 2.8 km s" 1 Mpc" 1 prior case tightens up the constraints by more than 1 a 
on J1a from below, while addition of SNIa data for the Hq ± er# = 73.8 ± 2.4 km s _1 Mpc -1 
prior case tightens up the constraints by more than 1 a on fl\ from above. Addition of SNIa 
data to the H(z) and BAO combination doesn't much improve the constraints on Q m o for 
either prior. 

Figures 191— TTTl show that adding SNIa data to the H(z) and BAO combination results in 
the most prominent effect for the XCDM case, Fig. [TOJ Here for the Hq ± cr# = 68 ± 2.8 km 
s _1 Mpc -1 prior it tightens up the constraints by more than 1 a on ux from above and below 
while for the H ± cx# = 73.8 ± 2.4 km s _1 Mpc -1 prior it tightens up the constraints by 
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Fig. 8.— Thick solid (dot-dashed) lines are 1, 2, and 3 o constraint contours for the 
0CDM model from a joint analysis of the BAO and SNIa (with systematic errors) data, with 
(without) the H(z) data. The full (empty) circle marks the best-fit point determined from 
the joint analysis with (without) the H(z) data (in the left panel the full and empty circles 
overlap). The a = horizontal axes correspond to spatially-flat ACDM models. In the left 
panel we use the H = 68 ± 2.8 km s _1 Mpc" 1 prior while the right panel is for the H = 
73.8 ± 2.4 km s _1 Mpc -1 case. For quantitative details see Table [3j 



more than 2 a on ujx from below. Addition of SNIa data to the H(z) and BAO combination 
doesn't much improve the constraints on Q m0 for either prior in this case. 

In the 0CDM case, Fig. [Til adding SNIa data to H(z) and BAO combination affects 
the constraint on a the most for the H ± o Hq = 68 ± 2.8 km s" 1 Mpc" 1 prior case. The 
effect on Q m Q is little stronger than what happens in the ACDM and XCDM cases but still 
less than 1 a. 

Table 4 lists the two standard deviation bounds on the individual cosmological param- 
eters, determined from their one- dimensional posterior probability distributions functions 
(which are derived by marginalizing the two-dimensional likelihood over the other cosmolog- 
ical parameter) for different combinations of data set. 



The constraints on the cosmological parameters that we derive from only the BAO an d 
SNIa data are restrictive, but less so than th ose shown in Fig. 4 of IChen fc Ratral (j2011bf ). 
This is because the new ISuzuki et all (120121 ) SNIa compilation data we use here is based 
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Model and prior 


v 2 . 
A mm 


B.F.P 


v 2 . 
A mm 


B.F.P 


v 2 . 
Amin 


B.F.P 


A 1 — \ TV X 

ACDM 
h = 0.68 ±0.028 


20.7 


QmO = 0.31 

Q A = 0.68 


566 


tt m0 = 0.30 
Q A = 0.70 


551 


fi m0 = 0.30 


/i = 0.738 ± 0.024 


21.0 


o n on 

Q A = 0.79 


567 


o n on 
\l m0 = 0.30 

Q A = 0.76 


tt A = 0.73 


XCDM 
h = 0.68 ±0.028 


20.7 


^mO = 0.31 
Ux = -0.99 


566 


Q-mO — 0.31 

Ux = -1.02 


551 


Q m0 = 0.30 


XCDM 
h = 0.738 ± 0.024 


20.8 


^mO = 0.28 
ux = -1-19 


567 


O m0 = 0.30 
co x = -1.08 


u x = -1.03 


0CDM 
h = 0.68 ± 0.028 


20.7 


QmO = 0.31 

a = 0.05 


566 


O m0 = 0.30 
a = 0.00 


551 


Q m0 = 0.30 


0CDM 
h = 0.738 ± 0.024 


22.0 


Q mQ = 0.29 
a = 0.00 


567 


Q mQ = 0.29 
a = 0.00 


a = 0.00 



Table 3: The minimum value of y 2 and the corresponding best fit points (B.F.P) which 
maximize the likelihood, for different combinations of data. The SNIa data values are for 
the case including systematic errors. 



on a more careful accounting of the systematic errors, which have increased. Consequently, 
including the H(z) data, in addition to the BAO and SNIa data, in the analysis, more 



signific antly tightens the constraints: compare Figs. [6]— 18] here to Figs. 4 — 6 of IChen fc Ratra 



( 12011bl ). We emphasize, however, that this effect is prominent only in some parts of the 
parameter spaces. 



7. Conclusion 

In summary, the results of a joint analysis of the H(z), BAO, and SNIa data are very 
consistent with the predictions of a spatially-flat cosmological model with energy budget 
dominated by a time- independent cosmological constant, the standard ACDM model. How- 
ever, the data are not yet good enough to strongly rule out slowly-evolving dark energy 
density. More, and better quality, data are needed to discriminate between constant and 
slowly-evolving dark energy density. 

It is probably quite significant that current H(z) data constraints are almost as re- 
strictive as those from SNIa data. The acquisition of H(z) data has been an interesting 
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Fig. 9. — Thick solid (dot-dashed) lines are 1, 2, and 3 a constraint contours for the ACDM 
model from a joint analysis of the BAO and H(z) data, with (without) the SNIa data. The 
full (empty) circle marks the best-fit point determined from the joint analysis with (without) 
the SNIa data. The dotted sloping line corresponds to spatially-flat ACDM models. In the 
left panel we use the Hq = 68 ± 2.8 km s _1 Mpc -1 prior while the right panel is for the Ho 
= 73.8 ± 2.4 km s" 1 Mpc" 1 case. For quantitative details see Table [3j 



backwater of cosmology for the last few years. We hope that our results will help promote 
more interest in this exciting area. Since the H(z) technique has not been as much studied 
as, say, the SNIa apparent magnitude technique, a little more effort in the H(z) area is likely 
to lead to very useful results. 
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